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Abstract 

The orbital Ml collective mode predicted for deformed clusters in a schematic 
model is studied in a self-consistent random-phase-approximation approach 
which fully exploits the shell structure of the clusters. The microscopic mech- 
anism of the excitation is clarified and the close correlation with E2 mode 
established. The study shows that the Ml strength of the mode is fragmented 
over a large energy interval. In spite of that, the fraction remaining at low 
energy, well below the overwhelming dipole plasmon resonance, is compara- 
ble to the strength predicted in the schematic model. The importance of this 

result in view of future experiments is stressed. 
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Among the collective excitations which may occur in metal clusters, the magnetic dipole 
mode predicted for deformed clusters in a schematic model [l] has unique and appealing 
properties which deserve a deeper investigation. This excitation, which is the analogue of 
the scissors mode predicted and observed 0] in deformed nuclei, is promoted by a ro- 
tational oscillations of the valence electrons against the jellium background. Indeed, in the 
semiclassical approach [[!]], the displacement field of the mode is composed of a rigid rota- 
tional velocity field plus a quadrupole term which comes from the boundary condition that 
the velocity flow vanishes on the deformed surface. The distortion of the momentum Fermi 
sphere generates a restoring force of the rotational oscillations. The mode is characterized 
by the magnetic quantum number K n = 1 + and falls at an excitation energy |l| 

1 20 7 

u m = V2u 5 7—— - -rN^'H eV (0.1) 
where ujq = {2epjmr 

2)1/2^-1/3 and 

u p are respectively the harmonic oscillator (HO) and 
the plasma frequencies, 5 is the deformation parameter, r s is the Wigner-Seitz radius, ep 
is the Fermi energy ( r s = 2.1 A and ep = S.leV for Na clusters), and N e is the number of 
valence electrons in a cluster. The latter is related to the number N of atoms in a cluster 
by N = N e or N = N e + 1 according that the cluster is neutral or has a positive charge 
Z = +1. The mode gets a Ml strength given by 

B(M1) = 2 | (K* = l + | X)k(z) | 0) | 2 /i' 
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where S = 2/3 N e m < r 2 > = 2/5 r 2 mN^ 3 is the collective mass parameter. As the formulas 
show, the Ml mode is peculiar of deformed clusters. Its occurrence would represent a unique 
and unambiguous fingerprint for the onset of quadrupole deformation. The main indicator 
of the deformation available so far is the splitting of the El resonance which, however, is 
often washed out or not properly resolved experimentally. The energy formula (1) reveals 
another appealing property. The Ml mode falls well below the energy of the overwhelming 
El resonance and, therefore, has good chances of being detected experimentally. In view 



of such a possibility, it is of the utmost importance to test the predictions of the schematic 
model by carrying out a microscopic calculation which fully exploits the shell structure of the 
clusters. Such a calculation should shed light on the microscopic mechanism which generates 
the mode and should reveal, eventually, new properties connected with the shell structure It 
certainly will ascertain if and to what extent the Ml strength is fragmented and quenched. 
Clearly such a mode can be detected and can be used as a signature for deformation only if 
its Ml strength remains concentrated in a reasonably narrow energy range. 

We performed our calculation in a self-consistent RPA (SRPA) approach which in 
the most general formulation |4j] is based on the Kohn-Sham functional [|7|. Here we skipped 
the self-consistent derivation of the one-body potential and adopted the phenomenological 
deformed Woods-Saxon well. On the other hand, we determined self-consistently the two- 
body potential starting with a set of displacement fields V/r,2i, where 



The resulting interaction was a sum of weighted separable terms peaked on different slices of 
the system. Due to its close link to the detailed structure of the system, such an interaction, 
in spite of its separable form, came out to be quite suitable for describing its dynamical 
properties. The SRPA fully exploits the shell structure and, as shown for the dipole response 
||, reaches the accuracy of the most refined and complete RPA approaches. Moreover, it 
preserves the simplicity of the schematic model |J which can be easily recovered if an 
anisotropic HO potential plus the single operator / 2 2i = r2 (^2i + ^21) are used. 

The choice of a quadrupole-like field was motivated by its close connection with the 
generator of the rotational oscillation, namely the angular momentum (see |J and Refs. 
therein). Such a link can be easily established for L = 2 in the HO space , where one finds 



being b(uj ) = 2u Q for the AAf = space and b(u ) = 5u for AJ\f = 2 (Af is a principle 
shell quantum number). On the other hand, the quadrupole fields act obviously also in 



fL2i = r L (Y 21 + Y 2 \), 



L = 2,4,6,8. 
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the E2 channel. Consistency requires that both, magnetic dipole and electric quadrupole, 
excitations should be treated contextually. It is worth noting that, the choice made for our 
displacement fields (see Eq.p.3|) is quite general for our purposes. Indeed, since spin-orbit 
coupling in clusters is negligible, the orbital excitations are decoupled from the spin ones, 
so that the spin-spin interaction can be safely neglected. The spin-quadrupole fields can 
be also ignored. In the nuclear systems, they are known to affect only the spin channel by 
renormalizing the spin-spin interaction. 

The parameters of the Woods-Saxon potential Vws — + exp(R(Q) — r)/a )] with 
R(Q) = Ro(Po + /?2^2o(@)) and Rq = r iV 1//3 were adjusted so as to reproduce the Kohn- 
Sham+SRPA results for the dipole plasmon in spherical sodium clusters |4j]. The fit yielded 
r = 2.5A, Vo = —7.2 eV and a$ = 1.25A for singly charged clusters and r = 2.4A, 
Vo = —5.7 eV and do = 1.11 A for neutral clusters. The values of the deformation parameter 
5 = ^4:5/1671^2 were extracted from the experimental data || (for N e < 34) following 
the prescription of Ref. ]]], or were taken from the calculations ]nj (for N e > 34). Only 
clusters with measured or predicted axial quadrupole deformation were considered. Equal 
deformation parameters were used for both charged and neutral clusters. As shown in Ref. 
||, the SRPA calculations with these parameters account well for the observed deformation 
splitting of the dipole plasmon in deformed Na clusters. 

The most meaningful results of the calculation are presented in Figs. 1-3. In Fig. 1 
the Ml strength distribution is plotted for singly charged clusters varying from N = 15 to 
295. In order to simulate the temperature broadening, we smoothed out the Ml strength 
with the Lorentz weight using the averaging parameter A = 0.05 eV. Such a simulation 
is in general rather rough as compared to an explicit treatment of electronic and ionic 
thermal fluctuations. Nonetheless, since we do not pretend to describe specific thermal 
effects, this simulation should be sufficient for our purposes, if we confine ourselves within 
the temperature interval 300-600 K, where jellium aproximation is appropriate. We also 
gave in Fig. 1 a quantitative estimate of the Landau damping by computing the width T 
of the resonance which ideally envelops all the peaks above a threshold value fixed to be 
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one-half the height of the highest peak ||. This definition yields the standard full width 
at half maximum (FWHM) in the simplest case of one-peak structure. The plot shows 
that, as the size of the cluster increases, the whole Ml strength is shifted downward with 
rising magnitude and fragmentation. Only in going from the light prolate Na^ to the light 
oblate Naf 5 , this trend is not observed. In heavy clusters with N e ~ 300, the Ml strength 
reaches the huge values 350-400 /i^. The fragmentation (Landau damping) gets also very 
pronounced, since T and u become comparable. Due to the small value of Co, however, the 
strength remains concentrated in a rather narrow energy interval. 

The softening of the mode as well as the enhancement of the Ml strength can be nicely 
explained within the semiclassical model with the decreasing importance of the surface with 
respect to the bulk as the sizes of the cluster increase. This causes a faster increase of 
the mass parameter with respect to the restoring force constant coming almost entirely 
from a surface shear, with consequent lowering of the energy centroid ( p.l|) and enhancement 
of the Ml strength ( p.2| ). A more detailed and exhaustive explanation is provided by the 
microscopic excitation mechanism. By expanding the deformed single-particle wave function 
into a spherical basis | m) = J2ni I nlm) and accounting for the fact that each | m) state 
is, in general, dominated by a single spherical configuration | nlm) , one obtains the transition 
amplitude 



Clearly, the main contribution to the transition amplitude comes from orbits with high 
angular momentum I and small magnetic quantum number m. On the other hand, orbits 
with high I values are present only in heavy clusters. Hence the enhancement of the Ml 
strength. At the same time, as the sizes of the cluster increase with consequent increment 
of the number of high values of I, the density of the particle-hole (p — h) levels increases and 
their relative spacings decrease, causing an overall downward shift and a more pronounced 
fragmentation of the Ml strength. The above formula enables to sharpen the geometrical 
picture of the mode. Since most of the strength comes from orbits with high I and small 
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m values, it follows that the oscillatory rotational motion is promoted mainly by the orbits 
which are almost orthogonal to the equatorial plane. 

In light clusters the Ml transition is promoted mainly by one or two configurations. 
Indeed, the left and right peaks are due by more than 95% to the p-h components [200]- [211] 
and [202]-[211] in Nat 5 , [312]-[321] and [310]-[321] in Nat 7 , and [321]-[310] and [321]-[312] 
in Na^ 5 , having adopted the Nilsson-Clemenger notation Nn z A |lTJ for the single-particle 



orbitals. The reason of the small p — h admixture induced by the residual interaction is 
simple. The p — h configurations are very few and far apart in energy. The only observable 
effect of the interaction is therefore a shift of the Ml strength. We may therefore conclude 
that in light nuclei the Ml mode has the character of a single-particle excitation. Only in 
heavy clusters the collective nature of the mode appears evident. 

Fig. 2 shows the Ml and E2 responses of Naf 19 over a much wider energy interval. For 
a more homogeneous comparison, we give the photoabsorption cross sections, a(Ml,fi = 
1) ~ ES(M,/i = l)u and cr(E2,fi = 1) ~ J2B(E2,/i = l)cu 3 , rather than the strengths. 
The Ml spectrum is composed of several, roughly equally spaced, resonances, enveloping 
closely packed transitions, coming respectively from AAf =0, 2, 4, ... p — h excitations. The 
group of AAf =0 transitions correspond to the low-lying Ml mode predicted in the schematic 
model [||]. The others have no classical counterpart. This large scale fragmentation limits 
drastically the extent of validity of approaches which rely entirely on sum rules. It is worth 
noting, on the other hand, that the cross section, being proportional to the energy weighted 
Ml strength, magnifies the high energy transitions. Had we plotted the Ml strength, we 
would have observed a most prominent peak positioned in the lowest energy region and 
several others, much less pronounced, at higher energy. It is also to be pointed out that 
the high energy peaks which are physically relevant, namely the ones below the ionization 
threshold (which is 3.8 eV in light clusters and 3.2-3.4 eV in the heavier ones), overlap 
mostly with the dipole plasmon resonance and, therefore, are hardly detectable. 

The lower panel of Fig. 2 shows that, consistently with the HO relation ( |0.4| ), the E2 
strength covers the same energy regions of the Ml strength. It is, however, dominant over 
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the Ml transition only in the intermediate region, which is in any case the domain of the 
dipole plasmon resonance, and almost absent in the low-energy region. This latter interval 
is exclusively covered by the Ml mode, consistently with the predictions of the schematic 
model. 

A more quantitative comparison with this model is presented for the low energy mode 
in Fig. 3. Although the schematic Eqs. (l)-(2) were derived for neutral clusters only, we 
considered both charged and neutral clusters. Indeed, the results change very little in moving 
from one kind to the other. The energy centroids of the low-energy Ml transitions scale 
with deformation and the number of valence electrons basically according to the law derived 
in the schematic model Jl|. The summed Ml strength scales according to the schematic law 
only in heavy clusters, but fluctuates strongly in the light ones. These fluctuations reflect the 
single-particle nature of the transitions and, in principle, invalidate the schematic model for 
clusters of these sizes. The low-lying Ml strength is of order of the semiclassical estimates 
and even larger in light clusters. Moreover, this strength summed in the interval 1-6 eV 
exceeds systematically (up to 10-100%) the estimate (2). This looks surprising since in the 
schematic model the mode gets the total Ml strengh by construction. This apparent paradox 
is solved if we recall that the semiclassical calculation of the strength is fully equivalent to its 
RPA evaluation in the AAf = HO space. Due to the degeneracy of the / configurations in 
the AjV = HO space, a given state | Nn z A) = J2i a v/ A I NIA) is not dominated by a single 
configuration | NIA), but involves contributions of all orbits with comparable amplitudes 
fljjfj . The resulting Ml transition amplitude does not get its main contribution from the 
orbits with largest angular momentum, as in our case, but is an algebraic weighted sum 
of different contributions (with small weights) from all configurations, with both large and 
small I values. Hence the enhanced Ml strength produced by our calculations. 

The main results of our RPA calculation are: i) The Ml strength, at least in heavy 
clusters, is not concentrated only at low energy, as predicted by the schematic model, but 
spreads over a large energy region among equally spaced peaks corresponding to AJ\f = 
0,2,4- • • p — h transitions, ii) In spite of that, closely packed Ml transitions still fall at 
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low energy and carry an overall strength which is comparable to the value predicted in the 
schematic collective model. Such a strength can become huge in heavy clusters. It can 
reach the impressive value of 350-400 fif already at iV e ~ 300. iii) The energy centroid 
of these transitions scales with the deformation and the number of valence electrons as 
in the schematic model only in heavy clusters. In the light ones, the strong fluctuations 
of the summed Ml strength invalidate the schematic model. iv)The crucial role of the 
quadrupole field in promoting the Ml mode is confirmed by the close correlation established 
quantitatively between Ml and E2 modes. While however the E2 strength is concentrated 
mostly in the uninteresting region covered by the plasmon dipole resonance, the Ml is the 
only dominant mode at low energy with a strength which becomes huge in heavy clusters. 
These properties render the mode quite accessible to experiments. Its occurrence not only 
would indicate the onset of deformation but would enable to measure the deformation itself 
by exploiting the scaling properties of the centroids with deformation and with the sizes of 
the clusters. 
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Figure captions. 

Figure 1. Energy distribution of the Ml strength over the interval 0-1 eV for sodium 
clusters ranging from N=15 to 295. The deformation parameter S, the energy centroid lD, the 
quantity T for estimating the Landau damping, and the summed Ml strength (J2B(Ml)) 
are given for each cluster. 

Figure 2. Plot of the Ml and E2 SRPA photoabsorption cross-sections over the full energy 
range in Naf 19 . The curves give the Lorentz averaged SRPA response, while the underlying 
bars show the pure discrete spectra, which better illustrate the Landau damping. 

Figure 3. Ratios between SRPA energy centroids (top) and Ml strengths (bottom), 
summed over 0-1 eV, and the corresponding schematic estimates (Eqs.1-2) for charged (stars) 
and neutral (triangles) clusters with N e =14,18,26,34,118,278, and 294. 
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